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ABSTRACT 

We cross-match galaxy cluster candidates selected via their Sunyaev-Zel’dovich effect (SZE) 
signatures in 129.1 deg^ of the South Pole Telescope 2500d SPT-SZ survey with optically 
identified clusters selected from the Dark Energy Survey (DES) science verification data. We 
identify 25 clusters between 0.1 < z < 0.8 in the union of the SPT-SZ and redMaPPer (RM) 
samples. RM is an optical cluster finding algorithm that also returns a richness estimate for 
each cluster. We model the richness d-mass relation with the following function (In dlMsoo) “c 
Bx lnM 5 oo + Ci ln£(z) and use SPT-SZ cluster masses and RM richnesses A to constrain the 
parameters. We find - 1.14^Qjg and Ci = Q.73^Qljy The associated scatter in mass at 
fixed richness is cri„^l^ = O.IS^^qqj at a characteristic richness A = 70. We demonstrate that 
our model provides an adequate description of the matched sample, showing that the fraction 
of SPT-SZ selected clusters with RM counterparts is consistent with expectations and that the 
fraction of RM selected clusters with SPT-SZ counterparts is in mild tension with expectation. 
We model the optical-SZE cluster positional offset distribution with the sum of two Gaussians, 
showing that it is consistent with a dominant, centrally peaked population and a sub-dominant 
population characterized by larger offsets. We also cross-match the RM catalog with SPT-SZ 
candidates below the official catalog threshold significance ^ = 4.5, using the RM catalog 
© 0000 RAS to provide optical confirmation and redshifts for additional low-^ SPT-SZ candidates. In this 

way, we identify 15 additional clusters with ^ e [4,4.5] over the redshift regime explored by 
RM in the overlapping region between DES science verification data and the SPT-SZ survey. 
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1 INTRODUCTION 


Clusters of galaxies were first identified as over-dense regions in the 
projected number counts of galaxies (e.g., |Abell|1958[[Zwicky et al.| 
[1968|l. Nowadays, clusters are also regularly identified through 
their X-ray emission (e.g., IGioia et al.|1990||Vikhlinin et al.|1998| 


Bdhringer et al.|2000 Pacaud et al.|2007'[ Suhada et al.|2012| > and 

at millimeter wavelengths through their Sunyaev-Zel’dovich effect 
(SZE) signatures ( [Sunyaev & Zerdovich|[T972) . Large, homoge¬ 
neously selected samples of clusters are useful for both cosmologi¬ 
cal and astrophysical studies, and such samples have recently begun 
to be produced using SZE selection ([Staniszewski et al.|2009|[Hi 


Iselfield et al.|20f^ [Planck Collaboration et al.|2015||Bleem et al. 


from optical and near infrared photometric surveys (e.g.. 

Gladders 

& Yee|2000 

IKoester et aL|20071|Eisenhardt et al.|2008l|N 

enanteau 

et al.|2010| 

Hao et al.|2010 Wen et al.|20121|Rykoff et al.|20141 

Bleem et al. 

2015a||Ascaso et al.|2014| and references therein), and 


even larger samples will soon be available from ongoing and fu¬ 
ture surveys like the Dark Energy Su rvey (DES,|The Dark Energy) 
Survey Collaboration 2005|‘ KiDS (de Jong et al.|2013|, Euclid 


I Laureijsrt^l j20lT i and LSST ( |LSST Dark Energy Science Col 
|laboration|2012| >. 

Reliable estimates of galaxy cluster masses play a key role 
in both cosmological and astrophysical cluster studies. Eirst, the 
abundance of galaxy clusters as a function of mass is a well-known 


cosmological probe dWhite et al.|19931|Bartlett & Silk 

19941 Eke 

et al.| 19981 

jViana & Liddle||1999[ Borgani et al.||2001 

Vikh inin 

et al.f20091 

Rozo et al.|2010||Mantz et al.|2010l|Allen et al.|201 11 

Benson et a 

.|2013 Boequet et al.|2015 and many others). Second, 


accurate estimates of cluster masses are crucial in disentangling 
environmental effects from the secular evolution processes shaping 
galaxy formation (|Mei et al.||20091 |Zenteno et al.||20111 |Muzzin| 
let al.|2012| >. 


In this paper, we calibrate the richness-mass relation for SZE- 
selected galaxy clusters detected in the DES science verification 
data (SVAl) using the redMaPPer (Rykoff et al. 2014) cluster¬ 
finding algorithm. Specifically, we study the clusters detected via 
their SZE signatures in the South Pole Telescope SPT-SZ cluster 
survey ( |Bleem et al.|2015'bl hereafter B15) that are also present in 
the redMaPPer catalog. We also study the distribution of offsets 
between the SZE derived centres and the associated optical cen¬ 
tres, properly including the SZE positional uncertainties. Einally, 
we demonstrate our ability to push to even lower candidate signifi¬ 
cance within the SPT-SZ candidate catalog by taking advantage of 
the contiguous, deep, multiband imaging available through DES. In 
this respect, our study points towards the combined use of DES and 
SPT datasets to provide highly reliable extended SZE-selected clus¬ 
ter samples. We note that historically the optical follow-up of SPT 
selected clusters was the original motivation for proposing DES. 

The plan of the paper is as follows. In Section]^ we describe 
the galaxy cluster catalogs and the matching metric we use in this 
work. Section describes the method we adopt to calibrate the 
SZE-mass and richness-mass relations. Our results are presented 
in Section]^ Sectioncontains a discussion of our findings and 
our conclusions. In the Appendix, we provide a preliminary analy¬ 
sis of a cluster sample created using an independent cluster finding 
algorithm — the Voronoi Tessellation (VT) cluster finder — which 
helps to highlight areas where the VT algorithm can be improved. 
Throughout this work, we adopt CIm = 0.3, Da = 0.7, Hq = 70 


' http://www.darkenergysurvey.org 


kms“* Mpc^', and erg = 0.8. Cluster masses are defined within 
R 500 , the radius within which the density is 500 times the criti¬ 
cal density of the Universe. Future analyses will include the de¬ 
pendence of the derived scaling relation parameters on the adopted 
cosmology by simultaneously fitting for cosmological and scaling- 
relation parameters (e.g. |Mantz et al.|2010[|Rozo et al.|2010[|Boc-| 
Iquet et al.|2015^ . 


2 CLUSTER SAMPLE 
2.1 SPT-SZ Cluster Catalog 

The SPT-SZ galaxy cluster sample used in this analysis has been 
selected via the cluster thermal SZE signatures in a point-source 
masked-region of 2365 deg’ of the 2540 deg^ (2500d) SPT-SZ sur¬ 
vey using 95 GHz and 150 GHz data. Typical instrumental noise 
is approximately 40 (18) /rKcMB-arcrnin and the beam EWHM is 
1.6 (1.2) aremin for the 95 (150) GHz maps. A multi-frequency 
matched filter is used to extract the cluster SZE signal in a manner 
designed to optimally measure the cluster signal given knowledge 
of the cluster profile, the SZE spectrum and the noise in the maps 
( |Haehnelt & Tegmark||1996[ |Melin et'^|2006^ . The cluster gas 
profiles are assumed to be described by a projected isothermal jj 
model ( jCavalie re & Fusco-Femian o|1976^ with P = \. Note that, 
as discussed in fVanderlinde et al.(^2010| ), the resulting SPT-SZ can¬ 
didate catalogs are not sensitive to this assumption. The adopted 
model provides a SZE temperature decrement that is maximum at 
the cluster centre and weakens with separation 6 from the cluster 
centre as: 

hT{e) = m[\ + {Bie,fr\ ( 1 ) 

where ALq is the central value and 6^ is the core radius. We adopt 
12 different cluster profiles linearly spaced from = 0.25 to 3 
aremin ( |Vanderlinde et al.|2010[|Reichardt et al.|2013| B15). For 
each cluster, the maximum signal-to-noise across the 12 filtered 
maps is denoted as f. The SPT-SZ cluster candidates with ^ > 4.5 
have been previously published in B15. 


2.2 DES Optical Cluster Catalogs 

The DES Science Verification Data (DES-SVAl) that overlap SPT 
have been used to produce optically selected catalogs of clusters. 
In Section [2.2.1| we describe the acquisition and preparation of the 
DES-SVAl data, and in Section [2.2.2| we describe the production 
of the redMaPPer cluster catalog used in the primary analysis. We 
remind the reader that in Appendix [A| we present results of a pre¬ 
liminary analysis of the VT cluster catalog. 


2.2.1 DES-SVAl Data 

The DES-SVAl data include imaging of ~ 300 deg^ over multi¬ 
ple disconnected fields jMelchior et al.|2015||Sanchez et al.|2014| 
IBanerji et al.|2015T >, most of which overlap with the SPT-SZ survey. 
The DES-SVAl data were acquired with the Dark Energy Camera 
( |Diehl T. et al.|201^ |Flaugher et al.|2012||2015| ) over 78 nights, 
starting in Fall 2012 and ending early in 2013 with depth compa¬ 
rable to the nominal depth of the full DES survey (Rykoff et al., in 
preparation). 

Data have been processed through the DES Data Manage¬ 
ment (DESDM, jPesai et al.||2012] > pipeline that is an advanced 
version of development versions described in several publications 
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( |Ngeow et al.|2006{ |Mohr et al.|2008| |2012| >. The data were cal¬ 
ibrated in several stages leading to a Gold catalog of DES-SVAl 
galaxies (Rykoff et al., in preparation). The Gold catalog covers 
~ 250 deg" and is optimized for extragalactic science. In particular 
it masks regions south of declination 6 = -61°, avoiding the Large 
Magellanic Cloud and its high stellar densities. Furthermore, the 
footprint is restricted to the regions where we have coverage in all 
four bands. 


2.2.2 redMaPPer Cluster Catalog 


The red-sequence Matched-Filter Probabilistic Percolation 
(redMaPPer, hereafter RM) algorithm is a cluster-finding algo¬ 
rithm based on the richness estimator of |Rykoff et al.|p012| l. RM 
has been applied to photometric data from the Eighth Data Release 
(DR 8 ) of the Sloan Digital Sky Survey (|Aihara et al.|201 H SDSS,) 
and to the SDSS Stripe 82 coadd data ( |Annis et~ |2014) , and has 
been shown to provide excellent photometric redshifts, richness 
estimates that tightly correlate with external mass proxies, and 
very good completeness and purity (|Rozo & Rykoff||2014 


|et al. 2014b|c^. We refer the reader to the paper by Rykoff et al. 


Rozo 


1 20 14^ for a detailed description of the algorithm. Here, we briefly 
summarize the most salient features. 

We employ an updated version of the algorithm (v6.3.3), with 
improvements summarized in |Rozo & Rykoff] p014[ l, |Rozo et al.| 
( |2015[ in preparation), and |Rykoff et al. j2015| in preparation). RM 
calibrates the colour of red-sequence galaxies using galaxy clus¬ 
ters with spectroscopic redshifts. RM uses this information to esti¬ 
mate the membership probability of every galaxy in the vicinity of 
a galaxy cluster. The richness A is thus defined as the sum of the 
membership probabilities (prm) over all galaxies: 


^ ^ Prm- (2) 

In addition to the estimate of membership probabilities, the RM 
centering algorithm is also probabilistic. The centering probabil¬ 
ity Pccn is a likelihood-based estimate of the probability that the 
galaxy under consideration is a central galaxy. The centering like¬ 
lihood includes the fact that the photometric redshift of the central 
galaxy must be consistent with the cluster redshift, that the central 
galaxy luminosity must be consistent with the expected luminos¬ 
ity of the central galaxy of a cluster of the observed richness, and 
that that the galaxy density on a 300 kpc scale consistent with the 
galaxy density of central galaxies. The centering probability fur¬ 
ther accounts for the fact that every cluster has one and only one 
central galaxy, properly accounting for the relevant combinatoric 
factors. These probabilities have been tested on SDSS DR 8 data 
using X-ray selected galaxy clusters, and have been shown to pro¬ 
duce cluster centres that are consistent with the X-ray centres ( |Rozo| 
|& Rykoff|20T4l ). 

The DES-SVAl RM catalog was produced by running on a 
smaller footprint than that for the full SVAl Gold sample. In par¬ 
ticular, we restrict the catalog to the regions where the z-band lOcr 
galaxy limiting magnitude is z > 22. In total, we use 148 deg" of 
DES-SVAl imaging, with 129.1 deg^ overlapping the SPT-SZ foot¬ 
print. In this area, the largest fraction (124.6 deg^) is included in 
the so called DES-SVAl SPT-E field. The final catalog used in this 
work consists of 9281 clusters with A > 5 and redshifts in the range 
0.1 < z < 0.9. Due to the varying depth of the DES-SVAl cata¬ 
log, RM produces a mask that determines the maximum redshift 
of the cluster search at any given location in the survey. As an ex¬ 
ample, the effective area in the SPT-E region at the highest redshift 
(z > 0.85) is only ~ 30deg^. In addition to the cluster catalog, the 


RM algorithm also uses the survey mask to produce a set of random 
points with the same richness and redshift distribution as the clus¬ 
ters in the catalog. The random points take into account the survey 
geometry and the physical extent of the clusters, and as with the 
clusters, only includes points that have < 20 % of the local region 
masked (see |Rykoff et al.|2014| >. 


2.3 Catalog Matching 

We cross-match the SPT-SZ catalog with the RM optical cluster 
catalog following the method of |Rozo et al.| ( |2014c^ . First, we sort 
the SPT-SZ clusters to produce a list with decreasing SZE observ¬ 
able and we sort the RM catalog to produce a list with decreasing 
richness. Second, we go down the SPT-SZ sorted list, associating 
each SPT cluster candidate with the richest RM cluster candidate 
whose centre lies within 1.5 R 500 of the SZE centre. Third, we re¬ 
move the associated RM cluster from the list of possible counter¬ 
parts when matching the remaining SPT-selected clusters. 

R 500 is first computed assuming the redshift of the optical 
counterpart and using the SZE-mass scaling relation parameters 
adopted in B15. We subsequently check that our sample does not 
change when adopting our best fitting scaling relation parameters 
(see Section [3TT] l. 

To test the robustness of our matching algorithm against 
chance associations, we first perform the above described proce¬ 
dure on a sample of randomly generated RM clusters as described 
in the previous Section. Positions of clusters in this randomly gen¬ 
erated sample do not correlate with the positions of the SPT-SZ 
clusters. Using an ensemble of 10'* random catalogs we measure 
the distribution of richness in chance associations for each SPT-SZ 
candidate. We then apply the algorithm to match the real RM clus¬ 
ter catalog with the SPT-SZ candidate list where f > 4.5. The distri¬ 
bution of probabilities p of chance associations estimated for each 
cluster candidate using the randomly generated samples is shown 
in Figure[^(filled histogram). 

In Figure (left panel) we show the resulting 84% and 95% 
confidence limits (solid and dotted lines, respectively) in the rich¬ 
ness distribution of the chance associations as a function of the 
SPT-SZ observable This test allows us to estimate the proba¬ 
bility of chance superposition for each SPT-SZ cluster candidate. 
As detailed below, we use this information to determine whether 
or not to include particular matches for further analysis. As Fig¬ 
ure]^ shows, this filtering of the matched sample then ensures that 
chance superpositions are playing no more than a minor role even 
at 4 < ^ < 4.5. 

Within the DES-SVAl region explored by RM there are 36 
such SPT-SZ cluster candidates. Using information on the contam¬ 
ination fraction at ^ > 4.5 of the SPT-SZ candidate list and on the 
redshift distribution of the confirmed cluster candidates ( |Song et ^ 
|2012| B15), we expect ~ 9 of these candidates to be noise fluctu¬ 
ations and ~ 80% of the real clusters to lie at z < 0.8. Therefore, 
assuming the optical catalog is complete, the expected number of 
real cluster matches is ~ 22. Similarly, one can estimate an ex¬ 
pected number of real cluster matches (23.6) by scaling the total 
number of confirmed clusters in B15 below z = 0.8 (433) in the 
2365 deg^ of the SPT-SZ survey by the DES-SVAl area overlap¬ 
ping SPT-SZ that has been processed with the RM cluster finder 
(129.1 deg^ ). 

The actual number of matches to the SPT-SZ candidates is 33. 
Eight systems are foreground, low-richness RM clusters that have 
been erroneously associated with SPT-SZ candidates with either 
previously measured redshifts (four systems) or lower limits esti- 
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Figure 1. Distribution of probability p of chance associations for the f > 
4.5 (filled) and f > 4 (empty) samples. Black lines show the sample used in 
this work, red lines show the sample rejected (Section[2.3). 


mated in B15 (the remaining four candidates) that are at z > 0.8; 
these systems are either noise fluctuations or real clusters that are 
at redshifts too high for them to be detected by RM. In fact, all of 
these systems have a probability p of chance associations estimated 
from the randomly generated sample that is p > 16%. Therefore, 
we remove these matches from the sample. This leaves 25 SPT-SZ 
candidates at ^ > 4.5 that have RM counterparts. We expect less 
than one false associations within this sample of 25 candidates. 
The associated optical richness as a function of the SPT-SZ sig¬ 
nificance is shown in the left panel of Figure This number is 
somewhat larger but statistically consistent with the expected num¬ 
ber of matches presented above. All 22 of the SPT-SZ confirmed 
clusters presented in B15 that lie at redshifts where they could 
be detected by RM are in this matched sample. According to our 
matching metric (which differs from the approach in B15), there are 
also three unconfirmed SPT-SZ candidates (i.e., candidates without 
identified optical counterparts in the B15 analysis) that have RM 
counterparts: SPT-CL 10502-6048, SPT-CL 10437-5307 and SPT- 
CL10500-4551. The newly confirmed clusters are highlighted with 
large circles in Figure]^ 

Of the 25 SPT-SZ candidates with robust RM counterparts, 
we use 19 of them to calibrate the RM richness-mass relation. 
Six clusters are excluded from the analysis for the following rea¬ 
sons. Clusters with estimated redshift z < 0.25 in the SPT-SZ 
catalog from B15 are highlighted in cyan in the left panel of Fig- 
ure[^ Because the j ^-mass relation is robust only above this redshift 
( [Vanderlinde et alT]|2010^ , these systems are not used in the fol¬ 
lowing analysis. Two systems (SPT-CL 10440-4744 and SPT-CL 
10441-4502) are excluded from this analysis as they are detected 
in SPT-SZ regions that have been masked due to their proximity to 
point sources, which can compromise the SZE signal-to-noise mea¬ 
surement. In addition, we exclude the three clusters highlighted in 
magenta: SPT-CL 10417-4748, SPT-CL 10456-5116 and SPT-CL 
10502-6048. These systems are strongly masked in the DES-SVAl 
data; based on the SZE position, the masks cover 40% of the to¬ 
tal cluster region. As a result, the associated optical counterparts 
are highly mis-centered, and the corresponding richness is severely 
biased. We note that the average centering failure rate caused by 
the detection mask is 12% (3 clusters out of 25), in comparison to 


the corresponding rate in the SDSS RM catalog, which is ==1% - 
2%. The difference reflects the fact that SDSS has a much larger 
contiguous area, while SVAl has a more aggressive star mask. We 
expect this failure rate will decrease as the DES coverage increases, 
and object masking improves. Furthermore, improvements will be 
made to the RM algorithm to estimate the masked area not only 
at the putative centre of the cluster, but at all possible centres. In 
this way, clusters at high risk of mask-induced mis-centering will 
be properly removed from the sample. 

The B15 catalog contains only SPT-SZ candidates with ^ ^ 
4.5. In this work we also apply the matching algorithm to SPT-SZ 
candidates at 4 < < 4.5. We identify 26 matches in this signal to 

noise range. The resulting probabilities of chance associations of 
the ^ > 4 sample are shown as empty histograms in Figure[T] Sim¬ 
ilarly to the f > 4.5 case, we exclude 11 of these systems, which 
have estimated probabilities p of chance associations p > 16%. 
For the 15 matched systems, the expected number of false associ¬ 
ations is also smaller than one. The remaining cleaned sample is 
shown as red points on the left panel of Figure]^ The resulting to¬ 
tal number of SPT-SZ and RM associations at ^ > 4 is 40. This 
number is in good agreement with the expectation (^ 36) obtained 
using the number of SPT-candidates above ^ > 4 in the DES-SVAl 
region explored by RM (88) and correcting it by the expected num¬ 
ber of noise fluctuations (~ 45) and the number of clusters above 
z > 0.8 (~7). We find that two ^ < 4.5 SPT-SZ candidates, SPT-CL 
10501-4717 and SPT-CL 10439-5611, have probabilities of random 
associations larger than 5%, and therefore it is not clear whether 
these low richness associations are correct (see Figure]^. 

The right panel of Figure contains a comparison of the 
redshifts from the RM catalogs with the redshifts published in 
B15 (zspt) for the same clusters (obtained through dedicated op- 
tical/NIR followup by the SPT team or taken from the literature). 
Clusters with spectroscopic redshifts are highlighted in green. We 
note that the redshift estimates are not biased for the clusters af¬ 
fected by masking (magenta points). For SPT-SZ candidates with 
4 < ^ < 4.5 the SPT collaboration did not complete followup opti¬ 
cal imaging, and therefore we adopt the redshifts of the RM optical 
cluster counterpart. 

Table [T] contains all SPT candidates with RM counterparts 
used in this work. For newly confirmed SPT-SZ clusters , the as¬ 
sociated Zspt redshift is not given. We caution that the masses for 
low-redshift clusters (z < 0.25) may be underestimated due to fil¬ 
tering that is done to remove the noise component associated with 
the primary CMB. 


3 MASS CALIBRATION METHOD 


We apply the method described in Bocquet et al. (2015) to char¬ 
acterize the T-mass relation of SPT-selected clusters. We refer the 
reader to the original paper for a detailed description of the method. 
A similar approach has been adopted by |Liu et^L]^2015) for study¬ 
ing the SZE properties of an X-ray selected cluster sample fr om 
the XMM-BCS survey i Suhada et al.|2012 Desai et al.j2012 1 . In 
this analysis, we consider the RM richness as a follow-up observ¬ 
able to the SZE-selected cluster sample. This choice is adequate as 


^ SPT-CL 10423-5506 was previously identified in Song et al. (2012) and 
Reichai'dt et al. (2012) at a redshift z = 0.21 ±0.04 with a signal-to-noise f = 
4.51. The associated redshift estimated with the same analysis presented in 
B15 is z = 0.25 ±0.036. 
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Figure 2. Left panel: Richness as a function of the SZE significance ^ for the matched cluster sample. SPT-SZ candidates with f < 4.5 (vertical line) are 
shown in red. Clusters at z < 0.25 (cyan) and clusters with miscentering due to a high masked fraction (magenta) are not used in the richness analysis. Large 
circles indicate the newly confirmed SPT-SZ candidates with f > 4.5. Solid and dashed lines represent the upper 84 and 95 percentiles in richness of chance 
associations of SPT-SZ candidates and clusters from the randomly generated RM catalog. Right panel: The estimated redshift for the RM sample as a function 
of SPT redshifts as presented in B15 from independent optical follow-up data. SPT-SZ candidates with spectroscopic redshift are shown in green. Magenta 
symbols are the same as in the left panel. 


there are no SPT-SZ candidates with ^ > 4.5 missing RM counter¬ 
parts in the redshift and spatial regime explored by the RM catalog, 
so that the cross-sample can indeed be thought of as solely SPT- 
selected. We note that this is not the case for SPT-SZ candidates 
with 4 < ^ < 4.5 that do not have RM counterparts. However, the 
adopted method is also accurate under the assumption that cross¬ 
matching the SPT-SZ candidate list with the RM cluster catalog 
cleans the SPT-SZ candidate list, removing the expected noise fluc¬ 
tuations. Within this context the resulting cluster sample is there¬ 
fore drawn from the halo mass function through the SPT-SZ selec¬ 
tion in the redshift range explored by the RM catalog. 

In the following subsections we describe the model we use to 
simultaneously constrain the SZE-mass relation (Section |3T} and 
the richness-mass relation (Section[3.2^. 


3.1 The SZE-mass Relation 

Following previous SPT papers (|Vanderlinde et aH2010[[Benson| 
|et al.|2013)|Reichardt et al.|2013[|fecquet et al.|2015| B15), we de¬ 

fine the unbiased SZE significance f as the average signal-to-noise 
a cluster would produce over many realizations of SPT data, if the 
cluster position and core radius were perfectly known. This quan¬ 
tity is related to the expectation value of ^ over many realizations 
of the SPT data by: 


^ = V<f>" - 3, (3) 

where the bias in (^) is due to maximizing the signal-to-noise over 
three variables (cluster right ascension, declination, and core ra¬ 
dius). The scatter of the actual observable ^ with respect to (^) is 
characterized by a Gaussian of unit width. The SPT observable- 
mass relation P(^\Msoo, z) is modeled as a log-normal distribution 


of mean 


<lnf|M5oo,z) 


InAszE + Rsze In 
E(z) 


M,, 


3 X 1014/7-1 Mf 


+CszE In 


E(z = 0.6) 


(4) 


and scatter Dsze, and where E(z) = H(z)lHo. At low significance 
^ < 2, there is a non-negligible chance of multiple low-mass clus¬ 
ters overlapping within the same resolution element of the SPT 
beam. We account for this by only considering the brightest of 
these objects per approximate resolution element and we compute 
following [Crawford et al.| ( |2010| l. The SPT observable- 
mass relation is therefore expanded to P(^max\OPi^\^soo, z) and 
(max is then converted to the observable ( as in Eq.3. 

To calibrate the (-M relation we use the subsample of clusters 
with ^ > 5 and z > 0.25 from the 2500d SPT-SZ catalog (B15). We 
determine the parameter values by abundance-matching the cata¬ 
log against our fixed reference cosmology. We predict the expected 
number of clusters as a function of mass and redshift using the halo 
mass function ( [Tinker et al.|2008) . We convolve this mass function 
with the observable-mass relation accounting for its associated un¬ 
certainties, and compare the prediction with the data. Our approach 
here is effectively the opposite of the typical analysis, where cos¬ 
mological parameters are deduced from the cluster sample using 
both priors and calibrating information to constrain the scaling re¬ 
lation parameters (e.g. [Benson et al.[[2013] [Bocquet et al.[[20T5l (; 
here, we assume perfect knowledge of cosmology to calibrate the 
scaling relation. Note that this method does not depend on any as¬ 
sumptions about hydrostatic equilibrium. 

We assume flat priors on Asze, Bsze, Csze and a Gaussian prior 
on £>sze = 0.18 ± 0.07; the latter corresponds to the posterior dis¬ 
tribution derived from the cosmological analysis of the full SPT 
sample (de Haan et al., in preparation). We obtain the following 
parameters for the f-mass relation by maximizing the likelihood of 
obtaining the observed sample in ( and redshift under the model 


© 0000 RAS, MNRAS 000, 000-000 































6 A. Saw, et al. 


Table 1. SPT-SZ cluster candidates with RM counteipaid. We report the SPT-ID (1), right ascension (2) and declination (3), SPT peak detection significance f 
(4), corresponding core radius (5), richness A (6), associated redshift from the RM catalog (7) and SPT catalog (8), and SPT derived masses (9). Coordinates 
are J2000. 


SPT ID 

R.A. 

DEC 


6c [arcmin] 

A 

ZRM 

ZSPT 

M5oo[10‘^/!7„'Mo] 

SPT-CLJ0438-5419 

69.574 

-54.319 

22.88 

0.50 

144.76 ±5.52 

0.42 ± 0.01 

0.42 

10.19 ± 1.33 

SPT-CL J0040-4407 

10.199 

-44.133 

19.34 

0.50 

137.45 ±7.03 

0.37 ± 0.01 

0.35 

9.71 ± 1.28 

SPT-CL J0417-4748 

64.344 

-47.812 

14.24 

0.25 

54.22 ± 6.75 

0.58 ± 0.01 

0.58 

7.41 ± 1.00 

SPT-CL 10516-5430 

79.149 

-54.510 

12.41 

1.50 

178.93 ±8.71 

0.33 ± 0.02 

0.29 

7.05 ± 0.97 

SPT-CL J0449-4901 

72.273 

-49.023 

8.91 

0.50 

91.37 ±4.75 

0.80 ±0.01 

0.79 

5.24 ± 0.78 

SPT-CLJ0456-5I16 

74.115 

-51.275 

8.58 

1.00 

73.08 ± 5.39 

0.56 ±0.01 

0.56 

5.39 ±0.81 

SPT-CL J0441-4855 

70.450 

-48.917 

8.56 

0.50 

86.96 ±4.55 

0.81 ±0.01 

0.79 ± 0.04 

5.10 ±0.77 

SPT-CL J0439-4600 

69.807 

-46.012 

8.28 

0.25 

55.18 ±3.52 

0.34 ± 0.02 

0.34 ± 0.04 

5.52 ±0.84 

SPT-CL J0440-4657 

70.229 

-46.964 

7.13 

1.25 

67.95 ± 3.62 

0.33 ± 0.02 

0.35 ± 0.04 

4.95 ±0.81 

SPT-CL J0447-5055 

71.843 

-50.921 

5.97 

0.25 

77.84 ±5.26 

0.40 ± 0.01 

0.39 ± 0.05 

4.24 ± 0.81 

SPT-CL J0422-5140 

65.591 

-51.674 

5.86 

1.00 

49.28 ±5.32 

0.58 ± 0.01 

0.59 ± 0.03 

3.98 ± 0.78 

SPT-CL J0439-5330 

69.928 

-53.502 

5.61 

0.75 

60.77 ±3.81 

0.43 ± 0.01 

0.43 ± 0.04 

3.97 ±0.81 

SPT-CL J0433-5630 

68.249 

-56.502 

5.32 

1.75 

60.75 ± 4.82 

0.71 ±0.02 

0.69 

3.56 ±0.78 

SPT-CL J0535-5956 

83.791 

-59.939 

5.20 

0.25 

50.25 ± 4.09 

0.67 ± 0.02 

0.58 ± 0.03 

3.46 ± 0.77 

SPT-CL J0440-4744 

70.242 

-47.736 

5.12 

1.25 

82.55 ±3.80 

0.30 ± 0.02 

- 

3.75 ± 0.83 

SPT-CL J0428-6049 

67.026 

-60.828 

5.11 

1.25 

55.91 ±5.95 

0.73 ± 0.02 

0.64 ± 0.03 

3.46 ± 0.79 

SPT-CL J0444-4352 

71.162 

-43.872 

5.01 

1.50 

70.53 ±5.91 

0.57 ± 0.01 

0.57 ± 0.03 

3.53 ± 0.82 

SPT-CL J0458-5741 

74.598 

-57.695 

4.87 

2.50 

37.90 ±2.68 

0.19 ±0.01 

0.19 ±0.02 

3.69 ± 0.85 

SPT-CL J0534-5937 

83.606 

-59.625 

4.74 

0.25 

40.43 ± 3.42 

0.58 ± 0.01 

0.58 

3.15 ±0.76 

SPT-CL J0502-6048 

75.724 

-60.810 

4.69 

0.25 

30.73 ±4.32 

0.79 ± 0.02 

- 

3.03 ± 0.76 

SPT-CL J0441-4502 

70.345 

-45.040 

4.62 

2.50 

51.22±4.11 

0.15 ±0.01 


3.49 ± 0.85 

SPT-CL J0429-5233 

67.430 

-52.559 

4.56 

0.75 

33.84 ±3.97 

0.52 ±0.01 

0.53 ± 0.03 

3.15 ±0.79 

SPT-CL J0452-4806 

73.002 

-48.108 

4.52 

0.50 

56.54 ± 4.89 

0.42 ± 0.01 

0.37 ± 0.04 

3.26 ±0.81 

SPT-CL J0437-5307 

69.259 

-53.119 

4.51 

0.25 

36.89 ±3.56 

0.29 ± 0.02 

- 

3.20 ± 0.80 

SPT-CL J0500-4551 

75.209 

-45.856 

4.51 

0.75 

32.68 ±4.82 

0.26 ±0.01 

- 

3.66 ±0.91 

SPT-CL J0453-5027 

73.307 

-50.451 

4.47 

0.25 

31.99 ±3.44 

0.77 ± 0.02 

- 

2.89 ± 0.74 

SPT-CL J0449-4440 

72.473 

-44.672 

4.37 

0.75 

54.50 ± 5.43 

0.15 ±0.00 

- 

3.42 ± 0.86 

SPT-CL 10423-5506^ 

65.809 

-55.104 

4.36 

1.25 

38.65 ± 3.44 

0.27 ± 0.02 

- 

3.26 ±0.83 

SPT-CL J0451-5057 

72.937 

-50.965 

4.34 

0.50 

83.11 ±4.62 

0.76 ± 0.01 

- 

2.81 ±0.74 

SPT-CL J0438-4629 

69.564 

-46.488 

4.31 

0.50 

41.01 ±3.37 

0.43 ± 0.01 

- 

3.07 ± 0.79 

SPT-CL J0456-4531 

74.099 

-45.523 

4.30 

0.25 

32.45 ±3.18 

0.29 ± 0.02 

- 

3.17 ±0.81 

SPT-CL J0431-5353 

67.970 

-53.896 

4.22 

0.50 

57.23 ±4.76 

0.75 ± 0.02 

- 

2.74 ± 0.73 

SPT-CL J0501-4717 

75.274 

-47.294 

4.20 

3.00 

19.47 ± 2.50 

0.35 ± 0.02 

- 

3.34 ±0.88 

SPT-CL J0518-5740 

79.507 

-57.670 

4.19 

0.25 

65.22 ±4.82 

0.82 ±0.01 

- 

2.60 ± 0.70 

SPT-CL J0438-4907 

69.655 

-49.117 

4.19 

1.75 

76.20 ±4.16 

0.24 ± 0.01 

- 

3.13 ±0.81 

SPT-CL J0513-5901 

78.273 

-59.029 

4.17 

0.25 

33.29 ±3.58 

0.61 ± 0.01 

- 

2.75 ± 0.73 

SPT-CL J0451-4910 

72.888 

-49.178 

4.14 

0.25 

54.10 ±4.10 

0.73 ± 0.02 

- 

2.71 ±0.73 

SPT-CL 10439-5611 

69.978 

-56.192 

4.14 

0.50 

17.82 ±2.93 

0.28 ± 0.02 

- 

3.10 ±0.81 

SPT-CL 10532-5752 

83.237 

-57.877 

4.11 

0.50 

48.49 ±3.81 

0.77 ± 0.02 

- 

2.59 ±0.71 

SPT-CL 10449-5908 

72.472 

-59.142 

4.11 

1.25 

107.14 ±5.29 

0.77 ± 0.01 

- 

2.68 ± 0.73 


derived from Eq. 3 and Eq. 4. The results are 


3.2 Richness-mass Relation 


AgzE “ 4.02 + 0.16, fiszE “ 1.71 + 0.09, 

CszE = 0.49 ± 0.16, Dsze = 0.20 ± 0.07. (5) 

For every cluster in the sample we also calculate the associated 
mass distribution, accounting for selection effects: 

z,f}oc P(^\M,oo,z,p) P(M,^\z, p), (6) 


As for the SZE-mass relation (Eq.|^, we assume a power law form 
for the T-mass relation: 


<lnT|M5oo,z) 


InA,! + Ba In 
-t-C,i In 


M500 


3 X 10i4/,-i Mf 
E{z) 


E(z = 0.6) 


(7) 


where the vector p encapsulates cosmological and scaling relation 
parameters and Pf^lMsoo, z,/7) is obtained from the ^-mass scal¬ 
ing relation as described by Eq. [^an dg The halo mass function 
P{MjQp\z, p) is the prior on the mass distribution at redshift z- 

Masses derived for the matched cluster sample are shown in 
Table 1. We note that both these masses and the SZE scaling rela¬ 
tion parameters quoted here are different from the ones reported in 
B15. We adopt the same fixed cosmology as in B15, but in this anal¬ 
ysis we consider data from the full SPT-SZ clust er-survey as op- 
posed to just the sample from the initial 720 deg“ 1 Reichardt et al. 
|20T3l >. 


where is the normalization, fi,i characterizes the mass depen¬ 
dence, and characterizes the redshift evolution. An additional 
parameter D,\ describes the intrinsic scatter in A, which is assumed 
to be log-normal and uncorrelated with the SZE scatter, with vari¬ 
ance given by: 

VardnTIMjoo) = exp(-<lnT|M5oo» + (8) 

The first term above represents the Poisson noise associated with 
the number of galaxies in a halo at fixed mass, and therefore we 
define the intrinsic scatter as log-normal scatter in addition to 
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Table 2. Best fitting parameters and 68% confidence level of the richness- 
mass scaling relation parameters described by Equation|7]and[^ 

Catalog __ Ih _O_ D.\ 

SPT-SZ+RM ^ > 4.5 66.1+“ 1.14+“i 0.73+“^ 

SPT-SZ+RM f > 4 69.8+“ 1-17+° 1-71+“ “ 0.20+“® 


Poisson noise. We assume flat priors on the distributions of A^, Ba, 
Ca and a positive flat prior for Da- 

The probability that a cluster with SPT-SZ signal-to-noise ^ is 
observed to have a richness A is 

PiA\^,z,ir) = jdM,ooPiA\M,oo,z,p)PiMsoo\^,z,ir). (9) 

The term P(A\M^qo,z, ^ contains the lognormal intrinsic scatter 
and normal measurement uncertainties in the observable A. We use 
the above distribution to evaluate the likelihood of the matched 
cluster sample defined through our cross-matching procedure. Note 
that we simultaneously vary both the optical and SZE scaling rela¬ 
tion parameters, further including the SZE data set from B15 with 
f > 5, z > 0.25 for constraining the SZE-mass relation. 
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4 RESULTS 

We present here the constraints on the richness-mass relation (Sec¬ 
tion |4T]( and then use these best fit parameters to explore whether 
the cumulative distribution of the matched samples are consistent 
with the expectations from the model (Section [4^ . Einally, in Sec- 
tion |4.3| we analysis the optical-SZE positional offset distribution. 

4.1 redMaPPer Richness-mass Relation 

We marginalize over the SZE-mass scaling relation parameters and 
constrain the posterior distributions for the RM T-mass scaling rela¬ 
tion. Our best fit parameters and 68% confidence level intervals are 
reported in Table and shown in Eigure|^ We note that the slope 
of the /1-mass relation is consistent with 1 within Icr (consistent, 
therefore, with the richness being proportional to the mass), and 
the model is consistent with no redshift evolution within Icr jAn-| 
|dreon & Congdon|2014^ . Eurthermore the resulting /1-mass relation 
is characterized by a remarkably low asymptotic intrinsic scatter, 
with cri„^|M 5 «, -» 0.15!““ as (In/flMjoo) oo (Eq.[^. Eollowing 
Evrard et al. (2014), we estimate the characteristic scatter in mass 
at fixed richness to arrive at (T\aM = 0.18!“* at /I = 70, 25% 

larger than the corresponding characteristic scatter in mass at fixed 
f. We present in Eigure|^the RM richness as a function of the SPT 
derived masses. Blue lines describe the best fitting model and in¬ 
trinsic scatter as derived from this analysis (Table 2) at a pivot point 
of z = 0.6. 

We have verified that our results are not dominated by uncer¬ 
tainties in the SZE-mass scaling relation by fixing these parameters 
to their best fit values. Our results are only marginally improved 
in this case. Consequently, future analyses with larger samples are 
expected to considerably reduce the uncertainties of the recovered 
/1-mass scaling relation parameters. 

In Eigure]^ one cluster appears to be an obvious outlier: SPT- 
CL J0516-5435 (^ = 12.4,A = 178.9,M5oo=7.05x10‘^/j“Mo). 
SPT-CL J0516-5345 is a well-known merger that is elongated in 
a north-south direction in the plane of the sky with an X-ray mass 
estimate nearly a factor of two larger than the SZE mass estimate. 
This cluster was in fact the strongest outlier in the sample of 14 


Figure 3. Posterior distribution for the four parameters of the /l-mass scal¬ 
ing relation (Equations |7] and [^. Predictions from the scaling relation of 
Rykoff et al. (2012) are shown as dashed black lines. Best fitting parame¬ 
ters and associated Icr uncertainties appear in Table|^ 

clusters in |Andersson et al!]j2011^ . |High et ar|p012| > made a weak 
lensing measurement of SPT-CL J0516-5345, and found that there 
was a significant offset between the brightest central galaxy (BCG) 
and the weak-lensing centre, consistent with the merger hypothe¬ 
sis. Additionally, [High et ah] p012| l found that the weak-lensing 
mass was in better agreement with the SZE-mass estimate than the 
X-ray mass estimate, at a level consistent with elongation observed 
in the plane of the sky. Therefore, SPT-CL J0516-5345 appears to 
be an outlier due to true intrinsic scatter in the observable-mass 
relations, so we leave it in our analysis. We note, however, that 
whether or not we include SPT-CL J0516-5345 in the fit has a sig¬ 
nificant impact on our results. Our best fit parameters shift from 
Ba - L14!“j‘andD.) = 0.1 5 !qj“ with this cluster, to fi,i = L00!“j 
and Da = 0.05!“3 when SPT-CL J0516-5345 is not included in 
the fit. Whether SPT-CL J0516-5345 represents a rare event in a 
non-Gaussian tail in the distribution of richness of galaxy clusters 
or the recovered log-normal scatter obtained when cluster SPT-CL 
J0516-5349 is included is more correct will thus need to await fu¬ 
ture analyses with larger samples. 

We convert the P(M^qq\A) scaling relation derived by [R ykoff] 
|et al.|j2012f using abundance matching and the SDSS maxBCG 
cluster-catalog ( [Koester et al.|2007] > to a richness-mass relation so 
that we can compare it to our results (see also |Evrard et al.|2014| (. 
The predictions from |Rykoff et al.H2012) are shown as dashed lines 
in Eigure under the assumption of no redshift evolution in the 
richness-mass relation. We note that all parameters of our derived 
RM /1-mass scaling relation for SPT-selected clusters are consistent 
with the |Rykofr et al.| ( |2012^ values. 

We repeat these analyses extending the sample to include 
those with 4 < ^ < 4.5 and find similar results (Table |^. The 
largest difference is in the redshift evolution term which now has a 
best fit value Ca = 1.71 !“*. While formally this difference does not 
have large statistical significance (L3 (t), it is coming from a sam¬ 
ple that includes a large fraction of the same clusters, so it is likely 
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Figure 4. Richness as a function of the SPT derived masses for the calibra¬ 
tion sample used in this analysis (Section 2.3). Blue lines show the best fit 
richness-mass relation and Icr intrinsic scatter. 


statistically significant. A larger redshift evolution term would im¬ 
ply that higher redshift RM clusters are less massive at fixed d. At 
the same time the derived scatter is also larger. 

Two of the clusters in the 4 < f < 4.5 range are compatible 
with false associations. Excluding the two matched clusters with 
the highest probability of random associations results in a ~ Icr 
shift in Bx (from Bx - 1.17 to = 1.04) and in a ~ 0.5cr shift 
in Cx (from Cx = 1.71 to Cj = 1.42), while the other parame¬ 
ters (Ax and Dx) are almost unchanged. Interestingly, even though 
we have increased the number of clusters in the SPT-SZ-l-RM sam¬ 
ple by 40%, the constraints on the scaling relation parameters are 
only mildly tighter. There are two reasons for this. First these lower 
signal-to-noise SPT-SZ clusters have larger fractional mass uncer¬ 
tainties in f ~ 0.23 and <^)“* ~ 0.13 respectively for the 

4 < ^ < 4.5 and ^ > 4.5 samples). Second, the richnesses are also 
systematically lower, leading to a larger Poisson variance. Thus, 
each low f cluster has less constraining power than a high ^ cluster, 
reducing the impact of extending the sample to include the lower 
mass systems. 


4.2 Consistency Test of Model 

We also test the consistency of the adopted scaling relation model 
with the data by examining whether we are finding the expected 
number of matches with the correct distribution in richness. To do 
this, we focus on the SPT-E field, which at ~ 124.6 deg^ is the 
largest contiguous region covered by the DES-SVAl data. We carry 
out two different tests. 

In the first test, we examine whether we are finding the ex¬ 
pected number of SZE-selected clusters and whether these clus¬ 
ters have the expected number of optical matches with the correct 
A distribution. We generate 10® Monte Carlo realizations of clus¬ 
ter samples extracted from the halo mass function above M 5 oo> 
lO'^ ® and assign richnesses A and SPT-SZ significance ^ us¬ 

ing the parameters we extract from our analysis of the real matched 
catalog. These Monte Carlo mocks are generated taking into ac¬ 
count the survey area as a function of redshift sampled by RM in 
the SPT-E field. We then apply the SZE selection-either ^ > 4.5 
or f > 4-and measure the cumulative distribution in A of the SZE- 


selected samples. We then compare this to the same distribution in 
the real matched catalog. 

Shaded regions in the upper panels of Figure]^ show 1, 2, and 
3cr confidence regions obtained from the mocks after marginalizing 
over the scaling relation parameters. The solid black line shows the 
distribution from the real catalog, which is in good agreement with 
the mocks. The largest observed difference is smaller than Icr indi¬ 
cating that the adopted model provides a consistent description of 
the observed number and richness distribution of the SZE-selected 
sample. Similarly, a Kolmogorov-Smirnov (KS) test for the ob¬ 
served cumulative distribution of matched systems as a function of 
A and the corresponding median distribution from the Monte Carlo 
simulations, shows that the null hypothesis of data being drawn 
from same distributions cannot be excluded and returns p-values 
0.76 and 0.96 for ^ > 4.5 and ^ > 4, respectively. 

The second test is focused on whether the RM cluster cata¬ 
log (which is significantly larger than the SPT-SZ catalog) has the 
expected number of SZE matches with the correct A distribution. 
Essentially, we take the observed RM catalog as a starting point, 
and calculate the expected number of systems with SPT-SZ coun¬ 
terparts given the model and parameter constraints from our T-mass 
likelihood analysis. This test differs from the first in that it takes the 
observed RM selected cluster sample as a starting point, and such 
a test should be more sensitive to, for example, contamination in 
the RM catalog. The formally correct way to evaluate a statistical 
difference between the RM selected sample and the SPT-SZ-l-RM 
matched sample used in this analysis would be to calibrate the A- 
mass and SZE-mass relations starting from the RM selected sam¬ 
ple. Such a study goes, however, beyond the scope of the current 
work and will be addressed in a future project. Therefore, the fol¬ 
lowing analysis is only intended to be a consistency check. 

For this purpose, we proceed by first computing, for each real 
RM selected cluster in the SPT-E field, the probability Pn, of that 
cluster also having ^ > 4.5 and therefore being in the matched 
sample. We define this probability as: 


Pm = 


pm,z)d^ 


- £ 

- 1 "^/ 


dMsoo P(^|A750o)P(4T50o|4, z). 


( 10 ) 


where P(M 5 oo|/l,z) oc PjAIMsoo, z)B(M 5 oo, z) and P(M 5 oo,z) is the 
halo mass function ( [Tinker et al.|2008) . 

To predict the expected number of RM clusters with SZE 
counterparts, we then randomly sample the scaling relation param¬ 
eters, determining for all the RM clusters in each case. We 
then Monte Carlo sample each P^ to produce randomly sampled 
matched cluster catalogs. We use the results from the ensemble of 
random matched catalogs to produce the expected cumulative dis¬ 
tributions in A. Orange, yellow, and red regions in the bottom panels 
of in Figure]^ show the 1, 2, and 3cr confidence regions, respec¬ 
tively, of the expected cumulative distribution in T for f > 4.5 (left) 
and f > 4 (right), given the RM catalog and scaling relation param¬ 
eter constraints as input. The dashed blue line shows the cumulative 
distribution of the entire sample of RM clusters in the area, while 
the solid black line shows the observed cumulative distribution of 
the real matched catalog. 

We note that the predicted number of SPT-SZ-l-RM matches in 
this case tends to be higher than that observed for A > 35, but the 
tension is weak. Of some concern is the high A end of the sample 
(A > 70), where only 9 of the 17 RM selected clusters have SPT-SZ 
counterparts at ^ > 4.5 despite their having large probabilities indi¬ 
cating they should be in the SZE-selected sample. However, a KS 
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Figure 5. Consistency tests of our model (Section]^ and the SPT-SZ+RM catalog from the ~124 deg^ SPT-E field. The solid black lines show the observed 
cumulative distribution in richness A of the f > 4.5 (left) and f > 4 (right) SZE-selected samples. Upper panels: Gray scale regions show 1, 2, and 3 it regions 
predicted by drawing 10® SZE-selected samples from the mass function and assigning A according to our scaling relation constraints. There is good agreement 
with the data. Lower panels: Dashed blue lines show the cumulative distribution in A of the full RM sample. Orange, yellow, and red areas define 1, 2, and 3cr 
regions representing the predicted cumulative distribution of the SPT-SZ+RM catalog using as input the full RM sample and the probability (Eq. |10| that each 
RM cluster will have an SPT-SZ counterpart, given our scaling relation constraints. 


test for the observed cumulative distribution of matched systems as 
a function of A and the corresponding median distribution from the 
mocks shows that the null hypothesis of data being drawn from the 
same parent distribution cannot be excluded and returns p-values 
0.90 and 0.33, respectively, for ^ > 4.5 and f > 4. 

While the KS test is showing no evidence for the two distribu¬ 
tions to differ, it is also known to be not very sensitive to the tails of 
the distributions jMoscovich-Eiger et al.|20l3) . With this respect, 
it is also therefore interesting to directly compare the observed cu¬ 
mulative distribution with the cumulative distribution from the pre¬ 
dicted number of matches. In this case, as for the classical KS test, 
we focus on the largest difference between the two distributions. 
Lower left panel of Figure]^ shows that at most, the observed cu¬ 
mulative distribution is in tension with expectations at the ~ 2cr 
level, providing some indication of tension between the observed 
sample and our model. For example, when restricting ourselves to 


clusters with A > 70, we find that our simulation results in a larger 
number of SPT+RM associations in 94.6% of our Monte Carlo re¬ 
alizations. This tendency to observe fewer matches than expected 
given the size of the RM selected sample could be explained by ei¬ 
ther contamination within the RM sample, additional incomplete¬ 
ness within the SPT-SZ sample beyond that caused by scatter in the 
SZF-mass relation, or simply a statistical fluctuation. Future work 
exploring the SZF properties of the lower mass systems along with 
an extension of the current analysis to the full overlap between DFS 
and the SPT-SZ survey will sharpen this test. 


4.3 Optical-SZE Positional Offset Distribution 

It has been shown that optical-miscentering can have a significant 
impact on the derived SZF signature for an optically-selected sam¬ 
ple ( [Biesiadzinski et al.|2012[[Sehgal et al.|2013]|Rozo & Rykoff| 


© 0000 RAS, MNRAS 000, 000-000 




















10 A. Saw, et al. 



Figure 6. SPT-CL J0433-5630: DES-SVAl gri pseudo colour image over¬ 
plotted with SPT-SZ signal-to-noise contours (in steps A^=l). The magenta 
circle shows the projected 7?5oo/2 radius at z = 0.69, while the green circle 
describes the Itr SPT positional uncertainty (Eq |l 1| . The cyan label marks 
the associated A ~ 60 RM cluster. 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Distance [R500] 


Figure 7. Solid histogram shows the measured fraction of SPT-SZ+RM 
clusters as a function of the optical-SZE positional offset in units of Rsoo- 
The green curve shows the SPT-SZ positional uncertainty, and the blue 
curves shows the best fitting SZE-optical positional offset model. 


[20T4l[Ro^et al.|2014a|bl l. We note, however, that in our case the 
SZE signal has been estimated at the SZE-determined position, so 
our results are not affected by optical miscentering. In fact, we can 
now use our data to constrain the distribution of offsets between the 
SZE-determined and the optically-determined cluster centres. 

As an example, we show in Eigure the DES-SVAl gri 
pseudo colour image of SPT-CL J0433-5630, an SPT-SZ selected 
cluster with ^ ~ 5.3 at redshift z = 0.69 (B15). Yellow contours 
show the SPT-SZ signal-to-noise in steps of A^ = 1, while the ma¬ 
genta circle describes the projected radius Rsool2. The cyan label 
refers to the associated RM cluster centre. This RM cluster has rich¬ 
ness A ~ 60. We note that the most probable central galaxy selected 
by RM is significantly offset from the SZE defined cenfre. As a re¬ 
sult, the measured SZE signature at the optical position = 4.1) 
would be significantly underestimated with respect to the derived 
unbiased quantity ^ = 5 obtained through Eq. 3. We stress that this 
effect is not important for the scaling relation results reported in 
Section [4T| as the sample analyzed here is SZE-selected. 

Eigure 1^ contains a normalized histogram of the distribution 
of cluster positional offsets in units of R 500 for the ^ > 4.5 analyzed 
SPT-SZ sample. Under the assumption that the measurement uncer¬ 
tainty from the optical side is negligible, we model this distribution 
as an underlying intrinsic positional offset distribution convolved 
with the SPT-SZ positional uncertainty. 

The 1(T SPT-SZ positional uncertainty for a cluster with a 
pressure profile given by a spherical p model with P = I and core 
radius d^, detected with SPT-SZ significance ^ is described by: 

= + ( 11 ) 

where ©beam = 1.2 aremin is the beam EWHM (see |Story et al.|20TT ] 
and |Song et al.|2012| for more details). As a result, the expected 
distribution of positional offsets in the case in which the intrinsic 
one is a (5-function is shown (arbitrarily rescaled) as a green line. 

|Song et ^ ( |2012) have shown that the intrinsic optical-SZE 
positional offset distribution for an SPT-SZ selected sample is con¬ 
sistent with the optical-X-ray positional offset distribution of X- 
ray selected clusters ( |Lin et al.|2004) . The offset distribution can be 
characterized by a large population of central galaxies with small 


Table 3. Best fitting parameters and 68% confidence level of the optical- 
SZE positional offset distribution. 


Catalog _po_ crolRsool crilRsool 

RM-^>4.5 0.63!°;>5 0.07+°;“ 0.25!«;O^ 


offsets from the SZE centres and a less populated tail of central 
galaxies with large offsets (e.g. |Lin et al.||2004| |Rozo & Rykoff| 
|2014[[Lauer et al.|2014| l. We therefore parametrize the distribution 
of positional offsets between the RM centre and the SZE centre for 
X as: 


P(x) - Ittx 


Po 

Inal 




1 -Po - 

2na\ 



( 12 ) 


where x = t/R^qq. While this model for the distribution was mo¬ 
tivated by the expected intrinsic positional offset distribution, the 
measured distribution will include both the actual physical SZE- 
central galaxy offset distribution and the systematics due to fail¬ 
ures in identifying the correct cluster center with the RM algo¬ 
rithm. Eor every cluster and parameter po 6 [0,1], (Tq 6 [0,1], 
and CT] 6 [(To, 1], we then compare the predicted offset distribu¬ 
tion obtained by convolving the model with the SPT-SZ positional 
uncertainty of Eq. 0 to extract the associated likelihood. Best fit 
parameters and 68 % confidence intervals are shown in Tablej^and 
joint and fully marginalized parameter constraints are shown in Eig- 
ure[ 8 ] 

We note that the positional offset distribution for the RM sam¬ 
ple is consistent with a concentrated dominant population (po = 
0 . 63 ^^ 25 ) of smaller offsets ((Tq = Rsoq) and a sub¬ 

dominant population characterized by larger offsets (ai = 0.25*111 
Rsoo)- In the limit where the SZ center and BCG are coincident 
for every cluster, the RM code( |Rykoff et al.|2012[ |2014^ predicts 
the SZ-RM offset distribution for this sample to contain a centrally 
peaked component with normalization (Pcen) to be 0.79. While this 
is formally consistant with the measured SZ - RM offset distri¬ 
bution, intrinsic optical-SZE positional offsets, such as those pre¬ 
sented by (e.g.,|Song et al.|[2012|) likely contribute to large sepa- 
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Figure 8. Posterior distribution for the 1 and 2cr level of the three parame¬ 
ter model describing the positional offset distribution of Equation^] Best 
fitting parameters are shown in Table[^ 


good agreement with those of |RykofF et al.|j2012| l and show that 
the scatter in mass at fixed richness A = 70 for this sample is only 
25% larger than the scatter in mass at fixed SPT-SZ observable 


2) Consistency test of model and matched catalog. We carry 
out two consistency tests to determine whether there is tension be¬ 
tween the observed matched sample and the expectations given the 
scaling relation model we have adopted. Both tests involve creat¬ 
ing Monte Carlo generated cluster catalogs with associated richness 
and SPT-SZ significance derived from the fitted scaling relations. 
The first test checks whether the correct number of SZE-selected 
clusters is found and whether those clusters exhibit the correct num¬ 
ber of optical matches with the expected A distribution. As is clear 
from Figure]^ the observations are perfectly consistent with the 
expectations from the model. Thus, our analysis shows that the 
data in our matched SPT-SZ-l-RM sample are well described by our 
adopted model. In the second test we take the much larger observed 
RM catalog as a starting point and use the model to test whether the 
expected number of SZE matches with the expected A distribution 
is found. Unlike the first test, this one would in principle be sensi¬ 
tive to contamination within the RM sample. Here the agreement 
is not as good because there is a tendency for there to be fewer 
observed matches than expected. However, the tension reaches the 
~ 2cr level at worst, and so there is no convincing evidence that our 
observed sample is inconsistent with the model. 


ration component of the observed distribution shown in Figure]^ 
Larger samples will be necessary to disentangle the impact of sys- 
tematics due to failures in identifying the correct cluster centre with 
the RM algorithm from the offset distribution due to cluster mor¬ 
phology. 


5 CONCLUSIONS 

In this paper, we cross-match SZE-selected cluster candidates with 
f > 4.5 from the 2500d SPT-SZ survey (B15) with the optically- 
selected cluster catalog extracted from the DES science verifica¬ 
tion data DES-SVAl. The optically-selected catalog is created us¬ 
ing the RM cluster-finding algorithm. We study the robustness of 
our matching algorithm by applying it to randomly generated RM 
catalogs. 

Using the adopted matching algorithm in the 129.1 deg“ of 
overlap between the two data sets, we create a matched catalog of 
33 clusters. Eight of these clusters are removed as likely chance su¬ 
perpositions that are identified using the randomly generated cata¬ 
logs. The resulting 25 cluster sample includes all previously known 
z < 0.8 and f > 4.5 SPT-SZ clusters in this area jSong et al.|2012| 
[Reichardt et al.|2013[ B15) in addition to three previously uncon¬ 
firmed SPT-SZ clusters. 

We then study three characteristics of this cross-matched SPT- 
SZ-hRM cluster sample: 

1) The richness mass relation of SPT-SZ selected clusters. We 
calibrate the T-mass relation from SZE measurements by applying 
the method described in Bocquet et al. (2015). In this analysis we 
assume a fixed fiducial cosmology and marginalize over the simul¬ 
taneously calibrated SPT-SZ ^-mass relation. We adopt flat priors 
on the richness-mass relation parameters. We find that the RM A- 
mass relation for SPT-SZ selected clusters is characterized by a 
small asymptotic intrinsic scatter = 0.15);qq® and by a slope 
= 1.14)^55'jg that is consistent with unity. Our constraints are in 


3) The SZE-optical positional offset distribution. We identify 
optical positional biases associated with 12% of the sample due 
to the masking in the DES-SVAl data. We remove these clusters 
and study the optical-SZE positional offset distribution for the rest 
of the matched sample. We model the underlying positional offset 
distribution as the sum of two Gaussians, while accounting for the 
SPT-SZ positional uncertainty. We show that the resulting distri¬ 
bution is consistent with being described by a dominant (63*2^%) 
centrally peaked distribution with (ctq = Psoo) and a sub¬ 

dominant (~ 37%) population characterized by larger separations 
(cTi = 0.25)^j5 gg Psoo)- For the same population, the RM algorithm 
assumes that 79% of the clusters will belong to a small-offset pop¬ 
ulation, consistent with our observations. 

We also match the SPT-SZ cluster candidates with 4 < ^ < 4.5 
to the RM optical cluster catalogs from DES-SVAl to extend the 
mass range of the SZE-selected clusters. Including the SPT-SZ can¬ 
didates between ^ = 4 and f = 4.5 increases the sample of matched 
clusters by ~ 40% compared to the f > 4.5 sample, highlighting the 
potential synergies of SPT and DES in producing lower-mass ex¬ 
tensions of SZE-selected cluster samples. We show that this larger 
sample produces results that are broadly consistent with the f > 4.5 
results, but only marginally tighter. This is due to the fact that mass 
constraints from lower signal-to-noise SPT clusters are somewhat 
weaker on a per cluster basis compared to the higher f sample. Fu¬ 
ture work benefiting from the larger region of overlap between the 
DES and SPT surveys will improve our derived constraints and help 
to better characterize the optical and SZE properties of cluster sam¬ 
ples in terms of positional offsets, purity, and completeness. More¬ 
over, the multiwavelength datasets available through DES and SPT 
enable characterization of the galaxy populations of large SZE- 
selected cluster samples, calibration of the SZE-selected cluster 
masses using weak lensing constraints, and many other promising 
studies. 
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APPENDIX A: STUDY OF MASS-RICHNESS SCATTER 
FOR CLUSTERS FOUND WITH THE VT METHOD 

The scatter of the mass-richness relation is a primary source of sys¬ 
tematic uncertainties in cosmological measurements using galaxy 
clusters. To ensure that this and other astrophysical systematics are 
kept under control in future cosmology analyses, the DES collabo¬ 
ration has pursued the development of multiple cluster finding al¬ 
gorithms. Here we present initial results obtained using the analysis 
framework described throughout this paper on a DES-SVAl cluster 
catalog created using the Voronoi Tessellation (VT) cluster finder 
( [Soares-Santos et al.|2011| l. The VT method is fundamentally dif¬ 
ferent from RM. Specifically, VT uses photometric redshifts to de¬ 
tect clusters in 2 - 1-1 dimensions, and is designed to produce a cluster 
catalog up to z ~ 1 and down to mass M ~ 10*^'^ Mq without any 
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assumptions about the colours of galaxies in cluster environments. 
VT has been tested on DES simulations ( |Soares-Santos et al.|201 1^ 
and on SDSS data. The mean mass-richness relation has heen cali¬ 
brated using a stacked weak lensing analysis of the SDSS VT clus¬ 
ters ( [Wiesner et al.|20T5| l. In this appendix, we describe the first 
study of the scatter of the mass-richness relation using our analysis 
framework. This study demonstrates that VT, although not as ma¬ 
ture a cluster finder as RM, is a promising algorithm. Our report 
provides an assessment of the current VT performance and helps to 
identify areas where improvements should he made. Some of these 
improvements can he applied to other cluster finding algorithms, 
too. 


Al VT method 

To detect clusters with the VT method, we huild 2D tessellations 
in each photometric redshift shell and flag galaxies that lie in high- 
density cells as cluster members. The density threshold is set in a 
non-parametric way from the 2-point correlation function of that 
given shell. This takes advantage of the fact that the distribution of 
VT cell densities can be uniquely predicted for any given point pro¬ 
cess. The 2-point function is a good description of the point process 
of the background galaxies on the sky. Clusters cause a small devi¬ 
ation from the predicted distribution, and we take the point where 
that deviation is maximized as the threshold for detection. 


A2 VT catalog for DES SVAl data 

For the DES SVAl data, the final catalog consists of 12948 clusters 
with richness N,,, > 5 and redshifts in the range 0.15 < z < 1. N„ 
is defined as the number of member galaxies. The catalog covers 
the SPT-E and SPT-W regions of the SVAl total footprint. We use 
DESDM data products as inputs, namely the Gold galaxy catalog, 
plus photometric redshifts and mask information. We used a mask 
to apply magnitude cuts 10 < mag_autoJ < 23.5. The photometric 
redshift information was obtained using a neural network method 
( [Sanchez et al.|20T4) . 

We match the VT and SPT-SZ catalogs using the same method 


as described in Section 2.3. We sort the SPT-SZ cluster sample ac¬ 
cording to decreasing SPT observable ^ and sort the VT catalogs 
according to decreasing richness. Then we associate the SPT clus¬ 
ter candidate with the richest cluster candidate whose centre lies 
within 1.5 R 500 of the SPT-SZ centre. We finally remove the asso¬ 
ciated optical cluster from the list of possible counterparts when 
matching the remaining SPT-selected clusters. This procedure re¬ 
sults in 42 VT clusters matched to ^ > 4.5 SPT-SZ clusters. 

A3 Results 

In Figure [AT| we show the optical richness Nyt as a function of ^ 
(left panel) for the 42 matched clusters. The scatter in the richness- 
f plot indicates that the VT richness performs poorly as a mass 
indicator. Improvements to the method are being developed based 
on these findings. Specifically, future work will explore using the 
galaxy magnitudes to calculate total stellar masses, which can then 
be used as a mass proxy. 

Figure 1^ also shows the VT estimated redshifts versus the 
redshifts determined in SPT follow-up observations (right panel). 
There is good agreement for z > 0.35. Deviations at lower z are 
found to arise from problems in the calculation of the 2-point func¬ 
tion predictions at z < 0.35, and a fix is underway. 

We also obtained the best fit parameters, and the correspond¬ 
ing 68% c.l. uncertainties, for the richness-mass scaling relation 
described by Equations |^and[^ 

Avt =48.1!“, SvT = 0.56!°i‘, 

CvT =-0.51!];*2, Dvt = 0.64!°;®. (Al) 

This result is consistent with the calibration for the mean relation 
obtained with weak lensing ( jWiesner et al.|20T5) (. The uncertain¬ 
ties, however, are larger than those obtained for RM. We expect that 
improvements to the richness estimator will result in better perfor¬ 
mance in future applications of the VT method. VT is not as mature 
a cluster finder method as RM but the complementarity of the two 
techniques argues for further development of this alternative. This 
study allowed us to identify a key area for improvement and es¬ 
tablishes a framework for future assessment of the mass-richness 
scatter for VT clusters. 
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Figure Al. Left panel: Richness as a function of the SPT-SZ significance ^ for the SPT-SZ+VT cluster sample. Clusters at z < 0.25 are not not used for the 
richness-mass fit and are shown in cyan. Right panel: The estimated redshift for the VT sample as a function of SPT-SZ redshifts as presented in B15. SPT-SZ 
candidates with spectroscopic redshift are shown in red. 
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